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At the suggestion of Dr. C. I. Bliss and by the courtesy of Professor 
H. G. O. Holck, whose data I shall use, the following note may serve 
to illustrate the extreme simplicity with which the technique derived 
from the analysis of covariance may be applied to problems concerned 
with the relation of a part to the whole, such as are constantly arising 
in many fields. 

The data consist of the body weights in kilograms and the heart 
weights in grams of 144 cats used in a group of digitalis assays: Of 
these 47 were females and 97 males. These data are presented in the 
table at the end of this note. To simplify the calculations only one 
decimal place was used for each value. Thus we have: 


TABLE 1 
ToTaL WEIGHTS 
Females Males 
Number 47 97 
Total body weight 110.9 Kg. 281.3 Kg, 
Total heart weight ............... 432.5 g. 1098.2 g. 
Heart as fraction of entire body -3900% .3904% 


The observed variation in these two measurements can, of course, be 
expressed by means of the sums of squares and products, as in the 
following tables. The rather intimidating phrase ‘‘spurious correla- 
lation’’ used in the earlier literature sometimes prevents workers from 
taking the simplest course. Obviously it would be easy to derive from 

1 Holeck, Harald G. O., Kazuo K. Kimura, and Barbara Bartels, “Effect of the 


Anesthetic and the Rate of Injection of Digitalis upon Its Lethal Dose in Cats,” Journal 
of the American Pharmaceutical Assn. 35: 366-370 (1946). 


the crude figures the corresponding square and product for the differ- 
ence between our variates, representing the weight of the body less the 
weight of the heart. 


TABLE 2 
SUMS OF SQUARES AND PRoDUCTS FoR BopDY AND HEART WEIGHTS 

Females df. (Body)2 (Body: Heart) (Heart )2 
ARG) HT Ras cste Santeannnte at aor 47 265.13 1029.62 4064.71 
Correction for Mean 1 261.677 1020.516 3979.920 
Dafkerences gecesi 46 3.453 9.104 84.790 

Males 
bal pre oy tne ete 97 836.75 3275.55 13056.17 
Correction for Mean . . 1 815.77 3185.07 12435.700 
iD sr siferdes SEN pean eae terme meee mes Nae 96 20.98 90.48 620.470 


From the corrected sums of squares and products we may find the 
regressions of heart weight on body weight, namely .2637% for females 
and .4313% for males. It will be noticed that these values are in the 
first case less and in the second case more than the average contribution 
of the heart to the total weight. The significance of such differences 
is often of importance, and must often appear to present a rather com- 
plex problem. 

We may, however, recognize in the above table the requisite data 
for the simplest form of an analysis of covariance, and test at once 
the homogeneity of the regressions therein. We may set up the re- 
gression equation H=a+bB, where H represents the heart weight 
and B the total body weight. If 6 is actually less than the average 
percentage, .89%, then a must be significantly greater than 0 and vice 
versa for b greater than .39%. a will be significant if the residual 
variance is significantly reduced by the inclusion of a in the regression 
equation. The partition of the variation in heart weight after adjust- 
ment for body weight by covariance furnishes such a test. Thus for 
females, since (1029.62)?/265.13 = 3998.481, we have a remainder of 
66.229 g.? for 46 degrees of freedom. Consequently, the significance 
of the observed difference between the regression value for heart 
weight and the average value can be tested by the following table. 


TABLE 3 
VARIATION IN HEART WEIGHT ADJUSTED FOR BoDY WEIGHT 
Females Males 
af. (Heart)2 M.S. af. (Heurt)2 MS, 
From original totals .00000... 46 66.229... : 96°) 233.669 7 
After correction for mean .......... 45 60.787 1.351 95 230.259 2.424 
Reduction due to mean oo... 1 5.442 5.442 1 3.410 3.410 
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The reduction due to the mean for females is almost significant (the 
variance ratio is 4.03 while the 5% significance level is 4.06), indicating 
a strong probability that a is actually different from 0. Hence in 


TABLE 4 


ORIGINAL DATA ON LivE BoDy WEIGHT IN KILOGRAMS AND FRESH HEART 
WEIGHT IN GRAMS OF FEMALE AND MALE DomEsTIC CATS 
USED IN EXPERIMENTS ON DIGITALIS 


Females 

Body Heart Body Heart Body Heart Body Heart 

Weight Weight Weight Weight Weight Weight Weight Weight 
2.3 9.6 2.0 7.4 2.1 9.8 2.9 10.1 
3.0 10.6 2.3 7.3 2.7 10.8 3.0 13.0 
2.9 919 2.2 7.1 2.2 9.1 2.2 8.7 
2.4 8.7 2.3 9.0 2.3 11.2 2.4 6.3 
2.3 10.1 2.1 7.6 2.1 8.1 2.4 8.8 
2.0 7.0 2.0 9.5 2.4 10.2 2.5 10.9 
2.2 11.0 2.9 10.1 2.7 8.5 2.5 9.0 
2.1 8.2 2.7 10.2 2.3 10.1 2.3 9:7 
2.3 9.0 2.6 10.1 2.1 8.7 2.3 8.4 
2.1 7.3 2.3 9.5 2.2 10.9 2.6 10.1 
2.1 8.5 2.6 8.7 2.3 7.9 2.3 10.6 
2.2 9.7 2.1 7.2 2.1 S35. ete. Peete; 

Males 

Body Heart Body Heart Body Heart Body Heart 

Weight Weight Weight Weight Weight Weight Weight Weight 
2.9 9.4 2.5 12.7 3.5 15.7 3.5 17.2 
2. 9.3 3.5 15.6 8 13.3 3.8 16.8 
2.2 7.2 2.4 9.1 2.2 yal 2.2 8.5 
2.9 11.3 2.2 7.6 2.5 7.9 3.3 15.4 
2.5 8.8 3.4 12.8 2.4 7.9 2.7 9.8 
3.1 9.9 2.6 8.3 3.9 14.4 3.2 11.9 
3.0 13.3 3.4 11.2 3.1 12:5 2.9 10.6 
2.5 12.7 2.6 9.4 3.7 11.0 3.6 13.3 
3.4 14.4 2.7 8.0 3.0 12.4 2.7 12.5 
3.0 10.0 3.3 14.9 3.2 13.5 2.9 11.8 
2.6 10.5 2.2 10.7 3.3 14.1 3.6 15.0 
2.5 8.6 3.2 13.6 3.0 12.7 2.8 10.2 
2.8 10.0 2.2 9.6 2.9 10.1 2.5 11.0 
3.1 12.1 3.5 11.7 3.0 10.4 2.6 11.5 
3.0 13.8 2.5 9.3 2.4 7.9 3.9 20.5 
2.7 12.0 3.2 12.3 3.8 14.8 3.0 12.2 
2.8 12.0 3.2 13.0 2.0 6.5 2.6 9.4 
2.1 10.1 2.7 9.6 3.1 11.5 2.7 9.0 
3.3 11.5 2.6 7.7 2.8 9.1 2.5 8.8 
3.4 12.2 2.7 9.6 2.3 9.6 2.2 9.6 
2.8 13.5 2.0 6.5 3.0 11.6 3.1 13.0 
2.7 10.4 3.1 14.3 2.2 (ee) 3.3 12.0 
3.2 11.6 2.4 7.3 3.4 12.4 2.7 itl 
3.0 10.6 3.6 14.8 3.5 12.9 3.6 11.8 
Rr Memerer A ce ny aii | Peder 2.8 11.4 


females heart weight probably increases less than proportionately to 
body weight. In the case of males, where the heart weight appears 
to increase more than proportionately to the body weight, the differ- 
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ence is not significant with these data. The close agreement between 
the sexes in the average percentage of the body taken up by the heart 
seems to mask a real difference in the heart weight to be expected for 
a given body weight. 

It may be noted that the estimated variance of heart weight for 
given body weight in males, 2.424 ¢.?, is considerably greater than the 
value for females, 1.851 g.2. The greater residual variance for males 
possibly was related to their larger size. The heaviest female weighed 
3.0 Ke. while nearly 40% of the males exceeded this weight. 
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2x2 FACTORIAL EXPERIMENTS IN INCOMPLETE GROUPS 
FOR USE IN BIOLOGICAL ASSAYS 


C. I. Buiss 


Connecticut Agricultural Experiment Station and Yale University 


Among the simpler experimental designs one of the most useful is 
the 2 x 2 factorial for testing two levels of two factors in all combina- 
tions. In developing an insecticidal treatment, for example, two con- 
centrations of an oil emulsion may be applied with two different 
stickers. In a biological assay one drug may be compared with 
another at each of two comparable dosage levels. The two-dose design 
has been used widely in assays based upon a graded response and pro- 
vides adequate estimates of relative potency when the form of the 
dosage-response curve is reasonably well known. 

The two-dose assay owes some of its popularity to the greater ease 
of assembling homogeneous groups of four than of six or eight re- 
sponses, which would be required for three or four dose assays. Even 
groups of four may be inconveniently large. When they consist of 
four successive reactions of the same animal separated by rest periods, 
the experiment is time-consuming and some individuals may die before 
it can be completed. When the groups are litters of the same sex the 
number with four males or four females may be limited. For these 
and other reasons groups of three or two may be preferred. Some 
experimenters [5, 6, 7] have used pairs and confounded one of the 
treatment effects with differences between pairs. This increases the 
error of the confounded effect. Frequently one would rather keep the 
same precision in all treatment comparisons and yet reduce the size 
of the group. 

The experimental design known as balanced incomplete blocks or 
groups meets this need. As originally described by Yates [9], only 
the information within groups was utilized. The method was adapted 
to two-dose assays arranged in pairs by Bliss and Rose [2]. However, 
groups with different treatment combinations also provide estimates of 
the effect of treatment. More recently Yates [10] has described sta- 
tistical methods for recovering inter-block information and has applied 
them to determining the adjusted mean response for each treatment 
and its error: In the 2x2 factorial design interest centers in three 
specific comparisons. Given this more limited objective, the calcula- 
tion can be reduced to a convenient form. 

Two restrictions in design. Frequently two restrictions are re- 
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quired in biological assays, such as occur in a Latin square. This is 
quite compatible with a balanced incomplete design. In a rabbit 
assay of insulin, for example, the responses of individual rabbits, the 
dates of testing and the doses of drug might form the rows, the columns 
and the letters of a 4x4 Latin square [1]. If the experiment were 
stopped after the third injection, the design would be an incomplete 
Latin square with the last column missing. Four rabbits would pro- 
vide a set, but additional sets could be added to gain the necessary pre- 
cision. In other cases the response may be modified by the number 
of times an animal has been used rather than by the environmental 
conditions on the day of testing, the columns then being the order of 
treatment, and the rows the individual animals. The 3x4 Latin 
rectangle, therefore, is a balanced incomplete design suitable for many 
biological assays. 

The analysis of incomplete Latin squares as described by Yates [8] 
is based entirely on the information within rows, columns and letters. 
It is thus subject to the same limitations as the original incomplete 
block design. The estimated effects of treatment in the two cases are 
identical and the error differs only by the segregation of the differences 
between the three complete columns. For this reason the calculation 
of the incomplete Latin square as originally described by Yates will be 
replaced by that in his later paper. 

The 2 x 2 factorial design in six balanced incomplete pairs can also 
be extended to segregate a second source of variation. Six pairs pro- 
vide all possible combinations of the four doses in an assay. If order 
of treatment is a potential source of error, the minimal number can 
be duplicated and its order reversed. Since even the basic set of 12 
pairs may not give sufficient replication, a second restriction on the 
design would seldom be burdensome. 

Thus two restrictions in design are compatible with 2x2 factorial 
experiments arranged in incomplete groups of three or of two units. 
Both numerical examples are of this type. A somewhat more elaborate 
table is needed for the original data and the analysis of variance has 
one more row than in assays having only the restrictions represented 
by the incomplete triplets or pairs. In other respects the analysis is 
the same. 

Experimental data. The method will be illustrated by two assays 
of parathyroid extract from the serum calcium of dogs [2]. The first 
example consisted originally of five 4x 4 Latin squares in which each 
row represented the responses of a single dog and each column one 


70 


of four different days of treatment. The letters designated two doses 
of 0.125 and 0.205 ce. per ke. of a ‘‘standard’’ preparation and two 
doses in the same ratio of an ‘‘unknown.’’ The data from the fourth 
test day have been omitted here, leaving the three observations on each 
of 20 dogs in Table 1. The rows in each Latin square have been rear- 


TABLE 1 


ASSAY OF PARATHYROID EXTRACT FROM 3 DETERMINATIONS ON EACH DoG OF THE 
SERUM CALCIUM IN MG. PER CENT (—10) 17 Hours AFrrrer INJECTION; 
DosaGce Ratio U,/U,=S8,/S,=1.640. Data or 
BLIss AND ROSE (2) 


Dose injected on Serum calcium, y 


Set 


3/15 3/24 4/5 3/16 3/25 4/6 Total, T, 


U, 8, 5.8 4.3 4.8 14.9 

Uz 8, 4.7 5.4 4.8 14.9 

Sy U, 5.0 4.5 4.0 13.5 

S, U, 3.8 7.0 6.0 16.8 

U;, 8, 7.0 6.5 5.0 18.5 

U, Sy 5.1 5.0 5.8 15.9 

S; U, 5.0 4.0 4.6 13.6 

S. Uz 2.0 3.8 4.0 9.8 

8, U, 3.6 5.3 7.2 16.1 

S, U, 4.4 3.8 4.4 12.6 

U, S, 5.8 5.0 5.2 16.0 

U, Sy 4.6 5.4 4.0 14.0 

U, Si 4.0 3.8 4.0 11.8 

U, S. 6.2 4.0 3.0 13.2 

S; U, 3.2 6.0 4.9 14.1 

8, Uz 3.0 4.0 4.0 11.0 

S, U, 3.0 3.4 3.8 10.2 

S, U, 5.8 6.0 5.0 16.8 

U, 8, 4.2 4.1 5.0 13.3 

U2 S, 5.2 6.2 5.0 16.4 

EEE) est A cee 914 97.5 94.5 283.4 

GLOUP occ a II Til IV Total, G 

Total, T, .... 71.5 73.4 70.5 68.0 283.4 


ranged so that the different combinations or groups of doses occur in a 
standard order. Since all dogs were treated at the same time, the 20 
observations on each day could be handled as a unit in segregating the 
variation between days. 

The second example required two injections in each dog separated 
by an interval of about two weeks. Twelve dogs comprised a basic 
set, the first six receiving the same treatment combinations as the sec- 
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ond six but in the reverse order (Table 6). Three such basic sets were 
completed, but unlike the Latin squares they were tested on different 
days with doses of 0.06 and 0.12 cc. of extract per kg. The dogs were 
more variable in size, from 6.4 to 16.4 kilograms. Except for propor- 
tioning the doses to the body weight, adjustments for differences in 
size have been omitted for the sake of simplicity. 


ASSAYS IN BALANCED GROUPS OF THREE 


Analysis of variance. This design is represented by the 3x4 
Latin rectangle with f replications of the basic set of 12 responses. 
An experiment may be replicated at different times, but here we will 
consider simultaneous treatments. Hach set is randomized indepen- 
dently by the procedure customary for Latin squares. We will adopt 
the convention that each row represents three responses by a single 
animal, each column one of the three different days of the test and 
each letter a different dose of drug. Two doses, S; and S2, are of a 
standard preparation and the other two, U, and U2, are equivalent 
doses of an unknown preparation, so far as its relative potency can 
be estimated in advance. The ratio of the two doses for the unknown 
must be the same as that for the standard. The designation of rows 
in any given case, of course, will depend upon the assay and often 
there may be no factor corresponding to columns. In all cases, how- 
ever, four rows or groups are needed for each basic set. 

It is convenient to arrange the rows of the basic set in a standard 
order for analysis. Each row will lack a different treatment and on 
this basis may be identified with one of four groups. The Roman 
numerals I to IV will designate groups without the treatments S2, 94, 
U, and U;, respectively. We then need the total for each column or 
day (7), for each row or animal within each set (7), for each of the 
four groups over all sets (7',), and for each dose separately within each 
group (Tz). The sums of each series of partial totals must check one 
another and equal the grand total (@). 

An analysis of variance is computed primarily to estimate two 
error variances, one within groups and the other between groups. The 
equations for computing the sums of squares are shown in Table 2. 
Those in the first four lines are based upon the marginal totals in the 
table of original data. From the individual entries we obtain the total 
sum of squares in line 8 and from the grand total the correction for the 
general mean (C), which is usually computed first of all. 

The above terms in the analysis of variance for our first example 
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TABLE 2 


CALCULATION OF THE SUMS oF SQUARES FOR THE ANALYSIS OF VARIANCE OF 
EXPERIMENTS WITH Two RESTRICTIONS WHEN ARRANGED IN f 
BAsiIc CONCURRENT SETS EAcH or Four TRIPLETS 


Mean square 


Variation due to aa is adres ae 

Columns (days) 8(T,7)/4f-C 

Individual rows (animals) - 8(T,*)/3-C 

Groups I to IV 3 8(T,?)/3f -—C 
Rows of same group .. 33 - Line 2—- Line 3 | 6,7+3 67 
Doses, ignoring groups .......... 3 8(V?) /12f 
Groups, eliminating dose ........ | 3 8(W?) /24f 6:7 + 80,7/3 
Intra-row error - Line 8—Lines | o? 

1, 4, 5 and 6 
fs 8(y*) -C 
Correction for mean .. Ra 1 C= 8?(y) /12f 
: : 4(3f-1) 
MbSE-TOW (OLTOR cosets ceca . Lines 4+6 o2 +———+ afr Oo? 


have been computed directly from the data in Table 1. The first step 
was to determine from the grand total the correction for the mean as 
C = 283.4?/60 = 1338.593. Substituting the partial totals in the for- 
mulae of Table 2 with f = 5, the sums of squares between days, between 
rows, between groups and the total followed easily. The results have 
been listed in Table 3. The sum of squares for groups I to IV ignor- 
ing dose was subtracted from that for rows to obtain the sum of squares 
between complete replicates or between dogs receiving the same com- 
binations of treatments. 


TABLE 3 
ANALYSIS OF VARIANCE FOR PARATHYROID ASSAY IN TABLE 2 


F Degrees Sum of Mean 
Variation due to i ae squares square 
Days (from T,) ... 0.930 0.465 
Rows (from T, 2 34.727 
Groups I to IV (from aie 1.011 
Dogs of same group (by difference) 33.716 2.1072 
Doses, ignoring groups (from V) .. 10.946 
Groups, eliminating dose (from W) # 1.941 0.6470 
Intra-row error (by difference) ........... 20.274 0.5793 = 8? 
Total (from y) ....... 67.807 
Correction for mean 1338.593 
Inter-row error 35.657 1.8767 =A’ 


The remaining two sums of squares are computed factorially. Fac- 
torial analysis may be characterized by factorial coefficients such as 
are given in the body of Table 4. The four treatments, U,, U2, S; and 
S, are shown at the top of the table and under each are three columns 
for the three groups (I through IV) in which it has been tested. 
There are thus 12 columns of coefficients. Below the coefficients is a 
row for the sums of the responses (7a) from the f replicates of each 
group within each treatment. Of the 11 degrees of freedom among 
these partial totals, six are isolated individually by factorial analysis. 
The remaining five form part of the intra-group error which is com- 
puted as a difference. 

The terms obtained factorially have been selected with two objec- 
tives in view, to complete the analysis of variance and to evaluate the 
assay. They are arranged in three pairs. The first line of each pair 
measures the crude effect of treatment, ignoring the subdivision into 
groups. The second line gives the net difference between groups 
associated with each specific treatment effect, after eliminating the 
effect of treatment. It is thus a component of the error. The treat- 
ments are the two direct effects and their interaction. Here they rep- 
resent (a) the comparison of the standard with the unknown, (b) the 
combined slope of the dosage-response curves for both standard and 
unknown and (ab) their interaction or the divergence of these two 
curves from parallelism. The factorial coefficients in each row are 
multiplied by the corresponding treatment totals (Z'z) in the last row 
of the table and the products summed. The sums of the products 
S(xTq) in the last column are designated as V and W respectively 
within each pair, with subscripts to indicate the treatment comparison. 

The factorial coefficients (2) meet two requirements. In each row 
their sum is zero and in any two rows the products of the correspond- 
ing coefficients add up to zero. Hence they are mutually orthogonal 
and the variances computed from the sums of products are indepen- 
dent of each other. In completing the analysis of variance in Table 2, 
the three values of V are squared, summed and divided by 12f ‘to 
obtain the sum of squares for doses ignoring groups. Similarly the 
three values of W are squared, summed and divided by 24f to obtain 
the sum of squares for groups eliminating dose. The values 12 and 
24 are the sums of the squares of the factorial coefficients in a given 
row or S(a?). The seventh line in Table 2 is the sum of squares 
remaining after subtracting those in lines 1, 4, 5 and 6 from the total 
in line 8. It measures the variation within rows and columns. The 
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inter-group variation in the tenth line totals the sums of squares in 
lines 4 and 6. 

Each treatment total (Tz) in Table 4 was obtained from Table 1 
as the sum of the five observations for a given dose in one of three 
groups. These totals were multiplied in turn by the coefficients of 
each line in the table and the products summed to obtain the listed 
values for V and for W. From the V’s in the first row of each treat- 
ment comparison—11.0, 23.0 and 2.6—we obtained the sum of squares 
(10.946) for doses ignoring groups and from the three W’s in the 
second row for each term the sum of squares (1.941) for groups 
eliminating dose. Hach represented three degrees of freedom. The 
analysis of variance was completed by computing the intra-group error 
as a difference and adding the two sums of squares for the inter-group 
variation. Each sum of squares was divided by the corresponding 
degrees of freedom to obtain the mean squares in Table 3. 

The two components of experimental error. The analysis of vari- 
ance in Table 2 provides estimates of two error terms, that within 
groups and that between groups. Our best estimate of the true error 
variance within groups (o,7) is the mean square, s,”, in line 7. The 
inter-group errors in lines 4 and 6 have been computed from two 
sources. The mean square in line 4 is based upon complete replicates 
and by component analysis [4] it is an estimate of o,?+380.?, where 
o2* is the population variance between groups. The mean square in 
line 6 is modified by the efficiency (H) of the experiment which here 
is equal to 8/9 as given by Yates [10]. Hence it is an estimate of 
o,7 + EH (802). 

The sums of squares in lines 4 and 6 are added in the last row of 
Table 2, giving a mean square designated as A’. Because the com- 
ponents from lines 4 and 6 have been weighted by their degrees of 
freedom, A’ is the composite in the last column of the table. Replac- 
ing the population variances o,? and o,? by their estimates, s,? and s.2, 


and solving for s,*, the net variance between rows can be computed 
from the mean square in line 10 as 


ep 4f—-1 

* 12f-4 
In experiments where the basic sets of 12 responses are not carried 
out together, the differences between sets should be segregated from 
the sum of squares in line 4. This would reduce the degrees of free- 
dom used as its weight in estimating s.2 and in consequence would 
modify Equation 1 and all terms depending upon it. 


EBay Od (EN ae ons (1) 
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If the arrangement in groups does not segregate a disturbing source 
of variation, the mean square between groups A’ should be equal to or 
smaller than that within groups (s,2). In this case treatment effects 
could be computed from the totals for each dose without considering 
differences between groups and the error variance would be sj’. 
Usually, however, A’ is larger than s,*._ Treatments are then adjusted 
for differences between groups and are subject to a somewhat larger 
error than if this adjustment were unnecessary. 

Numerical values have been obtained with the above procedure 
from the analysis of variance in Table 3. The mean square for intra- 
group error was s,?= 0.5793 with 35 degrees of freedom. The sums 
of squares for the variation among dogs were added and divided by 
16 +3=19 to obtain the mean square between groups of A’ = 1.8767, 
which exceeded s,2 significantly. Substituting in Equation 1 the esti- 


mated net variance between groups was re (1.8767 — 0.5793) = 


0.4402. The net variance between dogs was nearly as large as that for 
different readings on the same dog. 

Determination of adjusted effects of treatment. Two independent 
estimates of each treatment effect are obtainable from an experiment 
in incomplete balanced groups. One is based upon comparisons 
within groups and the other upon comparisons between groups. The 
most satisfactory measure of a treatment comparison is the weighted 
mean of the intra-group and inter-group estimates. The same result 
ean be computed alternatively and more conveniently from the values 
of V and W computed with the coefficients in any pair of rows in 
Table 4. These provide respectively a crude estimate of the treatment 
effect and a factor for its adjustment. The adjustment is first multi- 
plied by a coefficient », which is computed from the error variances, 
s2 and s.2. After a crude effect V is adjusted with » and W, the 
result must agree with the weighted mean from the intra-group and 
inter-group comparisons. 

The factorial coefficients in the first two lines of Table 5 lead respec- 
tively to the intra-group and inter-group treatment comparisons. Both 
represent the comparison of standard vs. unknown; the others follow 
asimilar pattern. The intra-group estimate (7';) is equivalent to that 
developed initially for the analysis of incomplete blocks [9]. Its 
variance may be computed as V(2;)'= fS(2?)s1 = 24fs,2. In units of 
a single response, the intra-group estimate is 


ti = T/fS |x| = T ;,/16f, 


which has a variance of 
Tif 


fS (x?) S51 
Vv ti = a = ° 
(4s) = Fegza] > B0F * a, 
The inter-group estimate, 7;, in the second line has an observed vari- 
ance of fS (a?) (s,?+ 852?) =12f(s,2+3s.2). In units of a single group 
of three responses, the inter-group estimate may be written as 


t; =T;/4f 


and its variance as 
V (ty )= 12f (s.? + 852?) /16f? = 3(s1? + 88,7) /4f = x. 
j 


These two estimates are averaged, weighting each by the reciprocal 
of its variance, to obtain the adjusted treatment effect 
wit; + W5tj 
te eereierna eR ME Ape (2) 
The terms ft; and t; in Equation 2 may be expressed in units of 
Vand W. The coefficients for the crude effect V and for the adjust- 
ment W have been copied from Table 4 to Table 5. The coefficients 
in the last two rows of Table 5 measure identically the same treatment 
effect as those in the first two rows. Thus the products of the corre- 
sponding coefficients in the first pair differ only in sign from those in 
the second. A second check is to compare the variances in the last 
column of the table, which were computed with the 7,’s of our numeri- 
cal example. The first two variances totalled 1.4741 + 0.6827 = 2.1568 
and the next two 2.0167 + 0.1401 = 2.1568. The factors for expressing 
7, and T; in terms of V and W ean be obtained by solving two pairs 
of simultaneous equations based upon the coefficients for U, and for S, 
in Group I. From the solution of these equations we find that 


T,=(4V+W)/3 and 7;= (V-2W) /3, 
which can be checked with the numerical values in Table 5. Substi- 
tuting in Equation 2 and simplifying we have 
is V 4 Wi Bw; Ww 
 12f wit w; 48f 


Each adjusted comparison tz, may be expressed conveniently in 
the same units as V by multiplying it by 12f. For the comparison 
of standard and unknown 


ta 


Taos fae (4) 
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Ww; — 8w; 
8(wit w;) 
factorial effect. In terms of the estimated variances within and be- 
tween groups, p= S2?/(3s,2+ 85.7). Substituting from Equation 1 


gue Efe ONG: eat) (5) 
asia (4f=1)A’+4(f —1)s,2° caused 


The three adjusted treatment comparisons may be determined for 
our numerical example. Substituting the values for s,? and A’ from 


c 19 x 1.2974 
> PB 152 x 1.8767 + 16 x 0.5793 


used in Equation 4 with the sums of products computed factorially in 
Table 4 to obtain the three adjusted effects of treatment of 7,=11.0- 
2 x 0.08370 x 4.1 = 10.314, 7, = 25.076 and Ta, =3.922. 

The error of the factorial comparisons. When the variation be- 
tween groups (A’) is larger than that within groups (s,”), the crude 
treatment comparisons are adjusted with the information between 
groups and the error for evaluating the adjusted comparison exceeds 
s,2._ In units of a single response the adjusted factorial comparison, 
t,, is the weighted mean of two components. Its variance is equal to 
the reciprocal of the sum of the weights or 


where p= and the subscripts a, b or ab refer to a given 


= 0.08370. This is 


Table 3 with f=5 


1 Sis iy =F 3S" 
V (ta) EF Wi Ss Wj > 4f : SCH an Osea z 
2 2 
Since 1+ p=), this may be reduced to 


V (ta) =i tae 


In practice, however, we are more interested in T, = 12ft, which has a 
standard error of 
Seen} (LES) sion eee (6) 


Each of the adjusted factorial comparisons, Tc, Ty, and Tap, has the 
same error. For testing their significance, a ‘‘just significant differ- 
ence’’ may be computed as 

C= Se Se ee eerie five rh eae (7) 


where ¢ is taken from a table of Student’s ¢ with the degrees of free- 

dom in s,?2 and the desired level of significance, usually at P = 0.05. 
Substituting the numerical values from our example, the standard 

error of each adjusted treatment effect was computed by Equation 6 
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as s = \/60 x 1.0887 x 0.5798 = 6.187. With 35 degrees of freedom the 
tabular value of ¢ at P=0.05 is 2.032, giving a ‘‘just significant differ- 
ence’? of c = 2.032 x 6.187 = 12.470 for evaluating the three adjusted 
factorial effects. In comparison with these criteria, the two dosage- 
effect curves were parallel (7.») well within the sampling error, the 
combined slope (7,) was highly significant and the potency of the 
unknown did not differ significantly from that of the standard (7,). 


THE BALANCED PAIR DESIGN 


Analysis of variance. Groups of two may be required in some 
experiments by the nature of the test material or by the exigencies of 
the assay. When each treatment is paired once with every other, six 
groups will cover all possible combinations of four treatments. How- 
ever, if the two responses in a pair represent successive reactions of 
the same animal, the six pairs of a series should be matched by a see- 
ond series in which the treatments are given in the reverse order. 
This leads to the experimental design shown in the left side of Table 6. 


TABLE 6 


ASSAY OF PARATHYROID EXTRACT FROM 2 DETERMINATIONS ON EACH DoG OF THE 
SERUM CALCIUM IN MG, PER CENT (—10) 17 Hourg A¥rrEer INJECTION; 
DosaceE Ratio U,/U,=8./S,=2.00. Data or BLISS AND ROSE (2) 


Serum calcium, y 


Treat- 
ment 


First set Second set Third set 


group | Ist 2nd 
T, |9/21 10/5) T, | 9/22 10/6 T, 
I 2.0 : 6.4 4.0 5.9 9.9 4.2 4.6 8.8 
Te 3.4 2.1 5.5 3.8 3.8 7.6 26 3.4 6.0 
2.2 2.2 4.4 2.8 4.0 6.8 34 4.8 8.2 
3.4 2.1 5.5 48 4.0 8.8 41 3.0 Tol 
24 4.2 6.6 3.4 3.6 7.0 3.0 4.0 7.0 
2.0 2.4 4.4 18 2.2 4.0 3.0) 76.2 8.9 
3.2 2.8 6.0 3.6 2.0 5.6 42 3.8 8.0 
1.0 2.5 3.5 Rat A) 5.8 3.0 3.3 6.3 
4.2 2.2 6.4 2.6 3.8 6.4 5.4 3.0 8.4 
16 2.4 4.0 1.6 4.0 5.6 40 8.8 2.8 
42 4.8 9.0 3.0 5.2 8.2 44 41 8.5 
3.4 1.6 5.0 5.2 4.0 9.2 46 3.2 7.8 
Total, T, and T, 33.0 33.7 | 66.7 | 39.5 45.4] 84.9 | 46.6 51.2 7.8 
OUD ce eten I Vit III DV ee VI Total, G 
Total, T, ... 44.7 34.7 40.6 43.8 46.3 39.3 | 249.4 


Even with 12 pairs in the basic set, f replicates may be required. 
While an assay in symmetrical pairs can be computed entirely from 
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| the differences within groups [2], the pair totals may furnish an 
appreciable amount of information. A more complete analysis differs 
but little from that for groups of three. In the following discussion 
we will assume that the basie set consists of 12 pairs and that there 
are f replicates of the basic set. 

The first stage in the calculation is to determine the column or 
day totals (T,), the row or pair totals (7,) and the totals for each of 
the six groups, (7). When the experiments are concurrent, the first 
four rows of the analysis of variance follow a pattern similar to that 
in Table 2. When the f basic sets of 12 pairs are not treated together, 
the sum of squares between complete replicates must be segregated as 
an additional term. It is computed from the f totals (7) of the dif- 
ferently-timed sets as shown in the first line of Table 7. Within basic 
sets the f differences between column totals (7’,) are independent of 
either treatments or rows. When the replicates are not concurrent, 
these totals are squared separately for each set, leading to a sum of 
squares between columns with f degrees of freedom (line 2, Table 7). 


TABLE 7 


CALCULATION OF THE SUMS OF SQUARES FOR THE ANALYSIS OF VARIANCE OF 
EXPERIMENTS WITH TWo RESTRICTIONS WHEN ARRANGED IN f BASIC 
Srtg of 12 Pamrg CarRiep Out aT DirrERENT TIMES 


De- 
A grees Mean 
Line els Sum of : 
Variation due to of square is 
No. free- squares estimate of 
dom 
1 | Complete sets or replicates | f—1 8 (7,7) /24-C 
2 | Columns within sets (days) | f 6(T,7)/12-C- 
Line 1 
3 Individual rows 1lf 8(T,*)/2-C- 
(animals) Ji rennees 5 Line 1 
4 Groups I to VI . e lif -5 8 (7,7) /4f - co 
5 | Rows of same group . e 3 Line 3 —- Line 677 + 20.7 
6 | Doses, ignoring groups ........ 3 8(V*) /24f 
Gand - 2 0, + 4027/3 
7 | Groups, eliminating dose ... 2 8[W?/2f8 (x?) ] 6.2 + 20,2 
8 | Untra-row CLVor. ceecccscssssscscees 11f-—3 | Line 9— Lines 
1, 2, 5, 6, 7 o? 
9 Total . Pe wine | 24f—-1 | &(y?) -C 
10 | Correction for mean . s 1 C= 8? (y) /24f pau 
TL | Tmter-row error ecscessssssssessesssen lif Lines 5+7 ort aes 0. 


Since it consisted of three basic sets tested at different times, the 
parathyroid assay in balanced pairs of Table 6 illustrates this second 
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type. The original values for each dog were summed to obtain the 
total for each pair (7',), for each day (7), for each set (7'.) and, over 
all three replicates, for each of the six groups (7,). The sums of 
squares based upon these totals were computed as shown in the work 
form of Table 7 with the results given in Table 8. 


TABLE 8 
ANALYSIS OF VARIANCE FOR PARATHYROID ASSAY IN TABLE 6 


Degrees 


hee Sum of Mean 
Variation due to ze “a squares square 
Complete sets (from 7's) 2 20.345. 10.1725 
Days within sets (from T,) . 3 2.353 0.7843 
Rows within sets (from 7,) a 33 45.301 
Groups I to VI (from Fg) nececenm 5 7.536 
Dogs of same group (by difference) 28 37.765 1.3488 
Doses, ignoring groups (from 7) ..... 3 27.183 
Groups, eliminating dose (from W) 5 1.317 0.2634 
Intra-row error (by difference) ......... 30 [ 18,403 0.6134 = s,* 
(otaln(iromity)) tancmenetoesar ts cea 71 107.366 
Correction for mean (from @) . sa il 863.894 
Enter-TOWserrons menue mae 33 39.082 1.1843 = A’ 


The sums of squares associated with the different treatments in a 
paired experiment are computed factorially with the coefficients in 
Table 9. The responses to each dose are summed separately for each 
of the six treatment combinations in groups I to VI, including the 
inverted pairs, and the totals (Zz) entered beneath the columns of 
factorial coefficients. The sums of squares in lines 6 and 7 of the 
analysis of Table 7 are computed with these coefficients. That in line 
6 is determined from the values of V and represents the crude effect of 
treatments or the differences in dose ignoring groups. As in the pre- 
ceding case, the coefficients (x) are assigned to the doses, Ui, U2, Sy 
and S, without regard to the pairs in which they occur. The three 
sums of products (V) are squared, summed and divided by 2f8 (2?) = 
24f to obtain the sum of squares for the analysis of variance. 

The sum of squares for groups eliminating dose has 5 degrees of 
freedom. Three of the five are associated with the three treatment 
effects ; their coefficients have been obtained from the equations for Ws 
as given by Yates [10]. The variation represented by the remaining 
two degrees of freedom is independent of that associated with treat- 
ment; their coefficients are shown in the last two rows of Table 9. The 
sum of products from each of these five rows (Wa to W-) is squared, 
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O'S VES SES SOL S9L B9T 63S GLB 6&3 PLT GBT PT ‘[e30} 980 
®=M=8E- ii T o- ii ga TL T Cin Cama 
PM=0'S T p= © T 0 fle i= il 0 (i 10119 SOp-10}UT 
Y= 8'8 - (i=" G qh T Coal i = 3 LS 
9 4=8'0 ta hee T T ea i T T i= aS preswg (q0) 
IM =9'L T a=! je ee he Z- 1 T (ure v3 
WA =FOF T ii T ts Se T if T (aS ~~ adotg (q) 
2M =9'8- 3 ie 4S 3 f= = T T 3 ii T 
PA=9L (Re bee jc SE eS T T T T T PE Tg GIN fl AGS) 
(Pr) The We san ae Ae Aca 1 ‘Nee 7 Nicos 0 joe Sint 
sponpoid *s *s “a ‘a WI9], 
yo wng 


dnor8 yove Ur syuswy}ve1} TOF (@) SJUOTOyFoos [VILOJOV MT 


} Wav, Wow vava ‘Z do Sdnowy NatMtag SHONGuaMETd Youd 
INANLSALTY ONV ININLVEYL, dO SLOTMIG AHL JO SISATVNY IVINOLOVT 
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divided by 2fS(x?) and the resulting variances added to obtain the 
sum of squares for groups eliminating dose. The remaining sum of 
squares, that for intra-pair error, is then obtained by difference from 
the total sum of squares, completing the analysis of variance. 

The mean square for intra-row error in Table 7 is an estimate 
(s:?) of the population variance o,2._ The two sums of squares which 
are totalled to obtain the mean square for inter-row error (A’) in the 
last line of Table 7 represent three components. The components from 
line 5 with (11f—5) degrees of freedom and from W, and W. (Table 
9) in line 7 with 2 degrees of freedom are based upon true replicates 
or their equivalent and estimate o,?+20,2. The third component, com- 
puted from W,, W, and Woz, is associated with the three treatment 
comparisons and hence is an estimate of o,?+ (20:2) where E = 4/6. 
The mean square, (A’), based upon the three components is a com- 
posite of the two error variances as shown in the last column and row 
of Table 7. Substituting the estimate s,? and s.? for o,” and op? and 
solving for s,*, the net error variance between rows is 


ALvny: 
“9 (Lif —1) 


The dose totals (Zz) from our numerical example in the last line 
of Table 9 were multiplied by the factorial coefficients in each row and 
the products summed to determine the required values of V and of W. 
The sum of squares for doses ignoring groups was calculated from the 
V’s as (12.6? + 42.4? + 0.8?) /72 = 27.183; that for groups eliminating 
dose was computed from the W’s in the remaining rows as (3.6? + 7.67 + 
8.8 + 3.87) /144 + 3.07/48 = 1.317. The intra-group error with 30 de- 
grees of freedom was obtained as a difference, leading to s,? = 0.6134. 
Finally, the two sums of squares for the variation between dogs were 
added to obtain the inter-group mean square in the last row of Table 
8, as A’=1.1843 with 33 degrees of freedom. The mean squares s,? 
and A’ differed significantly although less than in the preceding ex- 
periment. The error variance between rows was estimated by Equa- 


one as 3,3 = 3 (1.1843 — 0.6134) = 0.2944, 


82 (A’ — 5,2). eee 


Adjusted factorial effects and their error. The adjusted treatment 
effects from an experiment in balanced pairs resemble those for groups 
of three. Two independent estimates are available from the original 
data, one based upon intra-group and the other upon inter-group com- 
parisons. The first (¢;) is expressed in units of a single response and 
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the second (t;) in terms of a single pair of responses. By a process 
similar to that in Table 5, these estimates can be computed from the 
sums of products in Table 9 as 


ti = (2V + W)/48f and t; = (V - W) /12f. 
In determining the adjusted treatment effect t. by Equation 2, each 
component is weighted by the reciprocal of its variance or by 
w; = 16f/s,? and w;=2f/(s17 + 2s27). 
The adjusted effects ¢, are expressed more conveniently in the same 
units as V by the equation. 


sfV + 2 (sit + 4s*) W 9) 
9s? + 16s.” abet tet ae ereneres 


When the f basic sets are not concurrent, s.* may be replaced by its 


T= 24ft,=V + 


value in Equation 8 and the adjusted factorial effect computed as 


To=Viem{(L-k)Vt2(L tk) WY, nmin (10) 
where 
we Sit 28 11fA’—s,? be 
ie 95,2 + 16s,? 88fA’ + (lif 9)s at ae (11) 
and 
Qs,*___ lif(A*— 51’) (12) 


7 esr elif A’ se 5 


The standard error of 7, can be computed readily as 


ee [ean SN /O8Sh pein | eS As (13) 


which leads to c by Equation 7. 

The above equations apply directly to the numerical example in 
Table 6. The coefficients » and k were computed from A’ and s$,2 
: : a , 38.4686 18.8397 
in Table 8 with f=3, giving = 397 76g 7 OLN 1 and k = 35 7686 = 
0.48974. The adjusted treatment comparisons could be determined 
by substituting the values for V and W of Table 9 in Equation 10 or 


To. =Vo+0.1175 (0.5103V_ + 2.9795Wa). 


Then the adjusted values were T,=12.095, T,= 47.603 and I= 
— 2.233. By equation 13 each adjusted total had a standard error of 
s = \/864 x 0.11751 x 0.6134 = 7.892. With 30 degrees of freedom in 
s,2, the tabular value of t at P =0.05 was used in Equation 7 to obtain 
a ‘‘just significant difference’ of c = 2.042 x 7.892 = 16.115. When the 
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adjusted values were compared with these criteria, the experiment in 
balanced pairs led to the same conclusion as that in groups of three. 


THE ESTIMATION OF POTENCY 


The validity of a given assay may be judged from the relation be- 
tween the adjusted treatment effects and their error. The value for 
T, or slope should be highly significant and hence much larger than c. 
Ta», on the other hand, should be less than c and preferably less than s. 
This would indicate that the dosage-response curves for standard and 
unknown were parallel within the experimental error. A value of Tap 
larger than c and of T, equal to or less than s would suggest a qualita- 
tive difference between the standard and the unknown such that their 
relative potencies would vary with the level of response. However, if 
both 7, and Tq, were statistically significant, one might conclude more 
probably that the doses had exceeded the linear range of the dosage- 
response curve. In this case the experiment should be repeated with a 
revised assumption as to the potency of the unknown. As has already 
been noted, both examples met these requirements for a valid assay. 

The relative potency of the unknown may now be computed in the 
units which were assumed to be equal when setting up the assay. In 
terms of logarithms it is 

Meth /L a a | i een (14) 


where J, and J, are adjusted treatment effects and 7 is the interval 
between successive log-doses of drug. The proportionate relative 
potency or the anti-logarithm of M’ provides the basis for adjusting 
the concentration of an unknown preparation to equivalence with the 
standard. 

The standard error of M’, of use in estimating approximate confi- 
dence limits, may be computed as 


1 One a4 Y bs 
Sux eee ae Sa cel. eee (15) 

Approximate confidence limits are equal to 
Xe Ma toa ee ee eee (16) 


where ¢ is that used in computing c by Equation 7. With the small 
numbers in most assays, one would prefer (3) the exact confidence 
limits 


_ (ToT + VT? + Tr? - c?) 
~ Fa SRM a Bis AS 
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where all terms have their prior significance. The anti-logarithms of 
these limits give the range in the relative potency of the unknown at 
the confidence level selected in computing c by Equation 7. 

The relative potency of the unknown parathyroid extracts in our 
numerical examples has been computed with Equation 14. The ex- 
periment in groups of three used a dosage-ratio of 1.64, so that i = log 
(1.64) = 0.2148. The log-potenecy was determined as M’=0.2148 x 
10.3814/25.076 = 0.0883. The best estimate of the potency of the un- 
known, relative to the standard, was the antilog of 2+ M’ or 122.6 per 
cent. The equivalent values for the assay in balanced pairs with a 
greater dosage-ratio of 1=log(2) = 0.8010 were M’ = 0.0765 and rela- 
tive potency = 119.3 per cent. Hence the two unknowns agreed closely 
in potency. 

The approximate and exact confidence limits have been compared 
at P=0.05, using Equations 16 and 17 respectively. In percentage 
terms the approximate limits for the first example were 94 and 160 
per cent and the exact limits 96 and 180 per cent. In comparison with 
the exact values the approximate values for the first example under- 
estimated the limits of potency and the range. The two estimates 
agreed more closely in the second example, largely because the greater 
dosage interval increased the reliability of the slope. The approxi- 
mate and exact confidence limits were 94-152 and 94-158 per cent 
respectively. 

Summary. Factorial assays with two doses of the standard drug 
and two doses of the unknown have proved experimentally practicable 
and easy to analyze when based upon a graded response. In many 
cases potential sources of variation can be segregated by means of 
randomized block or Latin square designs. When homogeneous groups 
of four are less practicable than those of two or three, the balanced 
incomplete group may be substituted. The design and analysis of 
factorial experiments in balanced pairs and triplets are discussed and 
illustrated by two parathyroid assays based upon the serum calcium 


in dogs. 
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STATISTICAL TREATMENT OF COUNTS OF TRICHO- 
STRONGYLID EGGS* 


L. Ottis Emix** 
University of California, College of Agriculture, Davis 


Counts of nematodes’ eggs from fecal samples of the host have been 
used as a measure of the infestation level of the host although a multi- 
plicity of factors may affect the egg worm ratio. While techniques 
have been improved for securing reliable counts from samples, statisti- 
cal methods for analyzing these counts have seldom been used to full 
advantage. This is also true of other measures of infestation level. 

The ultimate aims of making estimates of infestation levels and 
statistical analyses of them vary with the experiment, but statistically 
speaking, the purpose is essentially the same, i.e., to show the presence 
or absence of significant differences between individuals, groups of 
individuals, or time intervals in the same individual or group of indi- 
viduals. The methods of analysis of variance could be applied to 
these types of data, in some cases without modification of experimental 
design. Gregory, Miller and Stewart [5], and Stewart, Miller and 
Douglas [7], tested differences between mean egg counts from sheep, 
breeds of sheep, different sires, and different ages by the individual 
pairing method, Snedecor [6] Section 2.13. Experimental design indi- 
cated the propriety of the method, but it was cumbersome and did not 
give the information afforded by analysis of variance. Ackert, Hisen- 
brandt, Wilmoth, Glading and Pratt [1], compared the resistance of 
five breeds of chickens to the nematode Ascaridia lineata by differences 
between means of number of worms and length of worms recovered at 
autopsy. Their data were suitable for analysis of variance and also 
covariance. More commonly, data on infestation levels have not been 
derived from experimental designs which would permit such analysis, 
but the experimental design usually could be altered to allow such 
analysis. 

The data presented here are typical, in most respects, of those which 
have been used for determining relative susceptibility of sheep to 
intestinal nematodes. This report presents a method of transforma- 
tion and an analysis of variance appropriate for the data. 

Twelve lambs were used to obtain these data. They ranged from 
157 days-of age to 180 days of age, with one 126 days old, at the end 


* Data taken from doctor’s thesis. 
** At present with the U. 8. Sheep Hxperiment Station and Western Sheep Breeding 
Laboratory, Dubois, Idaho. 


89 


of the experiment. Nine were Hampshires from two relatively new 
sires in the flock. The last three represented the Shropshire, Ram- 
bouillet and Merino breeds. 


TABLE 1 


ACTUAL COUNTS OF PARASITES’ Ea@cGs. MEANS FOR LAMBS, 
CORRECTED OR TRANSFORMED 


“sts Mean Corrected Mean square 
Bs eh ects square for root of 
root moisture corrected 
2510 8.3 2.73 12.9 3.36 
2563 5.5 2.22 7.6 2.62 
2566 6.1 2.27 ab i 3.08 
2569 13.9 3.61 19.8 — 4.33 
2571 15.8 3.82 19.7 4.31 
2574 17.9 4.16 27.5 5.11 
2578 11.5 3.33 ai 4.05 
2580 16.4 3.91 28.8 5.20 
2591 12.8 3.49 16.6 3.97 
2622 6.4 2.46 7.5 2.67 
2653 10.7 3.19 13.5 3.58 
BT76 8.9 2.89 10.6 3.12 
Mean 11.2 3.18 16.1 3.78 


These lambs were running with their mothers in the University 
flock at the beginning of the experiment. All Jambs were weaned and 
placed in one flock a few days later. Throughout the experiment the 
lambs grazed on Ladino clover during the evening and night, being 
corralled during the heat of the day. 

Counts of nematodes’ eggs were made over a period of 20 consecu- 
tive days. Fecal samples were collected each morning about 10 o’clock. 
The samples were prepared and counted according to the Stewart, 
Miller and Douglas [7] modification of the Caldwell method. Two 
counts were made on each sample. In making these counts a dilution 
technique was employed. One gram of feces was diluted to 10 ml. and 
two 0.1-ml. samples were completely counted. This is similar in prin- 
ciple to a hemacytometer count of blood cells, or the dilution method 
of counting bacteria. It is well established that such counts approach 
a Poisson distribution. 

Since these counts represented 12 heterogeneous individuals with 
considerable daily variation, one appropriate test of the fit of the basic 
data to a Poisson distribution was available, Fisher [4], p. 61. The 
two counts made on each sample were separate estimates of the same 
population of eggs, and a test of their variance would test agreement 


90 


with a Poisson distribution. Chi-square tests of the fit of individual 
counts to sample means should be distributed according to the proba- 
bility levels in the table of Chi-square, in this case with one degree of 
freedom. The Chi-square test, using the original uncorrected counts, 
Table 2, gave 6.25 with 7 degrees of freedom, a probability above 50 
percent. This test indicated that the data were distributed in a 
Poisson fashion. Further, the variances of the means were highly sig- 
nificantly correlated to the means. 


TABLE 2 
Trst or FIT or ORIGINAL UNCORRECTED CoUNTS TO A POISSON DISTRIBUTION 


Probabilit Observed Expected é 
limits value value Che sauate 
1.0-0.9 28 24 0.667 
0.9-0.8 20 24 0.667 
0.8—0.7 22 24 0.167 
0.7-0.5 38 48 2.083 
0.5-0.3 50 48 0.083 
0.3-0.2 23 24 0.041 
0.2-0.1 30 24 1.500 
0.1-0.0 29 24 1.041 
Total 240 240 6.250 
Probability of Chi-SqQuare o..ccccccscssmen 0.5 


The variability of such a series can be stabilized and the distribution 
made somewhat more normal by transformation to the square root of 
the count (actual number of eggs counted), Bartlett [2]. The mean 
square roots of the counts are listed in Table 1. The variances of the 
mean square roots were not significantly correlated to their means. 

Analysis of variance of the transformed variates was made accord- 
ing to Snedecor [6], Section 11.8. The method of making two counts 
on each daily sample afforded an estimate of sampling error. The 
sheep-day interaction, Table 3, was the appropriate term for experi- 
mental error. In this experiment, the sampling error was not appro- 
priate for use as the term for experimental error. Interaction was 
highly significant, denoting that oi.” was not 0, thereby increasing the 
variability between sheep beyond that due to sampling errors alone. 
Where greater time interva's with seasonal fluctuations were involved 
[5, 7] interaction was probably much less—judging from graphs of the 
data. The significant interaction indicated that the changes in the 
daily environment of the sheep or worms, or both, did not influence 
the concentration of eggs in the feces of all lambs in a similar manner. 
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The pairing method would, therefore, not have been appropriate for 
these data. 


TABLE 3 
ANALYSIS OF VARIANCE OF EGG CouNTS, TRANSFORMED TO THE SQUARE RooT 


Source Sum Degrees Estimate of variances 


Mean 
of of of panne 
variance squares freedom q Components Estimate 
Sheep siete 189.9033 ala! 17.2639" 6? + 20107 + 400,?_ —-o,? = .4093 
Day s.e ieteronatyeiceee 76.9052 19 4.0476* 67+ 201.7 + 240.2 oo? =.1316 
Sheep x Days ............ 186.0343 209 0.8901t 67+ 20157 O12" = .2796 
Within Subel. ....... 79.4504 240 0.3310 o o =.3310 


* Exceeds the .01 level using Sheep x Days as Exptl. Error. 
+ Exceeds the .01 level compared to Within Subclasses. 


The variances for the 12 lambs were found to be homogeneous by 
Bartlett’s test [6], Section 10.13. A standard difference, calculated 
from the error term, was thus appropriate for testing differences be- 
tween means. The standard difference between two means, at the 5 
percent level of probability was 0.585? eggs. Upon this basis, the 
mean counts of eight of the lambs were significantly higher than the 
lowest one. 

While the Chi-square test of agreement to Poisson distribution on 
the basic data showed agreement, Coehran’s test of the variance of 
the transformed variate showed non-agreement. Cochran [3] showed 
that 4x (within Subclasses sum of squares), using the square root 
transformation, should be distributed as Chi-square with Within Sub- 
classes degrees of freedom. For these data the test of significance, 
Fisher [4], Section 16, is \/2X?-—./2n—1 or \/2x4x79.4—\/2x240-1 
giving 3.3 which is highly significant. Even under conditions of this 
sort, Bartlett [2] has shown that the square root transformation is 
preferable if the variation of \/z is stable. 

Several known sources of sampling error remain in the technique, 
but cannot be eliminated. The ability of any worker accurately to 
discover all the eggs on every slide will not be perfect. In addition, 
the volume of sample placed on the slide will not be exactly 0.1 ml. for 
every slide. These were probably the main sources of error which 
raised the variance out of agreement with Poisson. 

Another source of increase in variance involves a correction for 
consistency of the feces, the factors being 1 for normally formed pel- 
lets, 1.5 for semi-mushy, 2 for mushy, and 4 for diarrhoeie feces. The 
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author’s data were originally in this adjusted form, Table 1, giving an 
estimate of co? of .518, and a Chi-square value of 9.7 with Cochran’s 
test. Since the adjustment for moisture content of the feces made no 
major changes in relative standings of the lambs, it would appear that, 
for these data at least, this correction should not be made. 

For many data on natural infestations, the additional correction 
for continuity, \/x+0.5, instead of \/x, would be necessary since the 
means are lower and counts of zero are not uncommon. While these 
transformations are time consuming, Bartlett [2] reminds the worker 
that ‘‘an unweighted average of the observed treatment response is not 
necessarily the best estimate of the true treatment response, and the 
average on the transformed scale will often be the better estimate when 
reconverted to the original scale.’’ 
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QUERIES 


(46) 


QUERY: In the formulas for missing values in experimental design, 
the error term is minimized. This of course inflates the / values and 
in turn the apparent significance of the effects and interactions under 
study. R.L. Anderson in the Biometrics Bulletin, Vol. 2, No. 3, gives 
formulas for removing this bias from the treatment-sum of squares. 
Would it not be safer in all cases to minimize the sum of squares due 
to effects regardless of what this did to the error mean square? 

This procedure has been used in one industrial experiment involv- 
ing four factors by simply placing ‘‘Y’’ equal to the lost value and 
equating to zero the first derivative with respect to Y of the total 
attributed sum of squares, four main ones, six first order interactions 
and four second order interactions. This appeared to give a satisfac- 
tory replacement for the lost value in this isolated case. 

What is the basis for preference in minimizing the residual mean 
square? 


ANSWER: Yates has shown in the articles referred to in my discus- 
sion of missing plot techniques that if a missing value is estimated by 
minimizing the sum of squares for error (SSH) the resulting error 
mean square is an unbiased estimate of o?. Even though the treat- 
ment mean square is slightly inflated, the bias is known and ean be 
removed if this is deemed advisable. It will be shown below that mini- 
mizing the sum of squares for treatments (SS7') will make SST too 
small (and in many cases far too small) and SSE too large. I want 
to present a short example of how misleading it would be to minimize 
SST for two treatments. Suppose we had these plot yields for two 
treatments : 
Replication 


Treatment Total 
1 2 3 4 5 6 
gL: 15 16 14 a7, 12 20 94 
2 10 9 9 13 8 y 49+y 
IDiMEreneS. caaascsnacsnenuks 5 7 5 4 4 20-y 454+y 


The minimum SST is obviously 0 and this is obtained by setting 
y=45. Hence we would have the first treatment 50% better than 
Treatment 2 in all five complete replications with the situation com- 
pletely reversed by the missing plot in replication 6. 
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The statement that by minimizing SST we would be safer ignores 
the facts that: 

(1) A test of significance should be as discriminating as possible, 
and a test based on a minimized SST (with resultant inflated SSH and 
far too small #7) could be quite unreliable, as it is in the example cited 
above. 

(2) Often the main use of the analysis of variance is to furnish a 
best unbiased estimate of o?. The method of estimating the missing 
value by minimizing SSH does not give an unbiased estimate of o”. 

Another point would plague one if he estimated y by minimizing 
SST. How could he handle the missing degree of freedom? Since 
SSE has been inflated, it does not seem reasonable to reduce the degrees 
of freedom for error, and it would be hard to explain that we had one 
less degree of freedom for treatments. This problem could become 
quite perplexing for several missing values. 

The foregoing discussion is clarified by considering the theoretical 
basis of the analysis of variance. The analysis of variance can be con- 
sidered as a special application of normal regression theory. I will 
illustrate the general principles for a simple randomized block experi- 
ment with ¢ treatments and r complete replications. The yield of any 
plot, yi;, can be represented as 

(1) Yig = B+ TE + pj t €i5 Ua etate los, freer La tagsetas 
where p is some fixed over-all effect, 7; is a fived effect of the 7-th treat- 
ment, p; is a fixed effect of the j-th replicate, and e;; is the residual 
effect, unexplained by the three fixed effects. It is assumed that ei; 
are normally and independently distributed with mean 0 and common 
variance o”. 

The best unbiased estimates of the fixed effects can be obtained by 
use of the method of least squares, whereby the residual sum of squares, 


(2) BRei;? = BE (yis— 7s — ps)”, 


is minimized with respect to each of these fixed effects, subject to the 
conditions, 


(3) Sse! 


Tf p, 7, pj are replaced by their best estimates, m, ¢; and rj, respec- 
tively, in equation (2), an unbiased estimate of the sum of squares for 
error is obtained : 


(4) SSE =33(yi;-m-—ti-7;)?. 
The mean square for error, s?, is an unbiased estimate of o°. 
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In order to test the null hypothesis that all the 7; are 0, the original 
estimating equation (1) is replaced by 
(Dy Yig= me tp it eaj 
and a new residual sum of squares SSE’ is determined by minimizing 
3Xe’:;? with respect to pw’ and the pj. The difference between SSH’ 
and SSE is an unbiased estimate of the sum of squares for treatments, 


SST, given by the analysis of variance. SSE and SST are indepen- 
dent estimates of ox?; hence, the null hypothesis can be tested by 


means of the Ff ratio, where F’ apd /S2: 


In the usual analysis of variance, SST is computed as 


37? G2 
(5) Een 


where 7; is the total yield over all replications for the i-th treatment 
and G is the total yield for the entire experiment. This SST is un- 
biased and independent of SSE only when all treatments are included 
in each replicate. If any treatment has a missing observation, an un- 
biased SST’, independent of SSE, can be obtained only by computing 
SSE’ — SSE, which have been determined for only the existing data. 
Any substitute method must be evaluated as to how closely it can pro- 
duce unbiased and independent estimates of SST and SSE. 

The missing-plot technique of estimating the value of a missing 
plot, say yu, by minimizing SSE furnishes an SST biased slightly 


upwards (by _ ‘t,?), However, we see that SSE is unbiased. The 
residual sum of squares is 
(6) 32’ (Yij —m —t,-13)? + (Yu —m—t,—14)2, 
where 33’ extends over only the existing data. When this residual is 
minimized with respect to yu, it is seen that y= m+ti+71; hence, 
the residual sum of squares is unchanged and is exactly SSE. 

If yi: were estimated by minimizing SST from equation (5), it can 
be shown that 


=m+t,+r ne 
Yu= 1 ae eh es 


Applying this result to (6), we find that in this case SSE is inflated by 


[ : AR Similarly, SST is too small by aus ree 
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For those who wish to use the method of covariance, I might men- 
tion that in this case estimating equations (1) and (1)’ merely have a 
new variable y,, added. The regression method then would adjust all 
sums of squares for y,:, so that unbiased and independent values of 
SST and SSE would be determined. 

It should be pointed out that these results could equally well be 
obtained by use of the method of maximum likelihood, and the general 
conclusions given here can be extended to more complicated experi- 


mental designs. 
R. L. ANDERSON 


(47) 

QUERY: The following problem has worried me for some time. In 
almost all modern textbooks of statistics the procedure given for com- 
paring two samples as to mean and variance involves the following: 

Consider two samples of NV, and N, variates respectively ; we wish 
to know whether they can have come from the same population. The 
squares of the deviations of the variates of each sample from its own 
mean are added together and divided by N1+N2-2. This is con- 
sidered an estimate of the variance (s?) of the population from which 
both samples were drawn, on the null hypothesis that they were in fact 
from a common population. Why is it not more accurate to obtain 
such an estimate of s? as follows: take the grand mean of all variates 
in both samples, take deviations from it, square, add, and divide by 
N,+N.-1? That would seem to be the obvious way to get a best 
estimate of combined variance, and usually it will lead to larger esti- 
mates than will the conventional method. 


ANSWER: Apparently you do not distinguish clearly two different 
problems, estimation and testing hypotheses. Perhaps also you have 
not discriminated among the circumstances in which each problem is 
to be solved. 

Consider first the estimation of variance. You have given the cor- 
rect method if the two samples are assumed to have been randomly 
drawn from a common normally distributed population; but I see no 
reason to think that in this case your estimate ‘‘usually will lead to 
larger estimates than will the conventional method.’’ This conven- 
tional method is appropriate if the sampled populations are assumed to 
have the same o? but not necessarily the same means. The correspond- 
ing analysis of variance is this: 
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Source of Degrees of Mean 


Variation Freedom Square 
al! My 
Individuals N,+N.-2 M, 
Total N,+N.-1 M, 


M, is here an estimate of o?; it is a weighted mean of the estimates 
given by the two samples, the assumption being that they are esti- 
mating a common population variance. My is now an estimate of 


P : Ni? +N,? ; 
a +hooy?, where ky is N,+N,-———., and where ox? is the com- 
Ni+N, 

ponent of variance introduced by the difference between the population 
means. Under these conditions your estimate, M7, does not estimate 
o* but will usually be larger than o? because of the component cy’. 
This is the situation you probably had in mind in making the last state- 
ment in your query. 

Let us now look at the conventional test of significance, F =My/M,; 
that is, the test of the null hypothesis, ox” = 0, under the assumption 
that o? is the variance common to the sampled populations. An essen- 
tial condition is that, under the null hypothesis, My and Mz shall be 
independent and shall both be estimates of o2. If you used Mp (your 
proposed estimate of variance) in the denominator, the resulting ratio 
would not be distributed as F’; moreover, it seems obvious that the test 
corresponding to your proposal would be less sensitive than F. 

For the t-test a similar argument can be constructed, the principles 


being the same. 
J. G. DArroca 


(48) 

QUERY: I have & samples of r with varying sample size. I wish to 
test the hypothesis that they are random samples from a common 
normal population. Using Fisher’s z-transformation, the quantity, 
(SNz)? 

Vile? 
is said to be distributed as chi-square with (k-1) degrees of freedom. 
T have tried to derive this by taking 


(1) SN2? — 


(2) x’ = 3% (0-e)*/e, 
but without success. Nor have I been able to derive the formula from 
(3) x? = (32?) /o?. 
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It seems to me that in (3) we have only the numerator of formula (1). 
That is, unless we can assume that 2 is normally distributed with unit 
standard deviation. Would you please put me on the right track? 


ANSWER: Formulas (1) and (3) are seen to be equivalent if the 
following pertinent facts are noted. (i) The quantity z is almost nor- 
mally distributed with variance 1/(n—3) where n is the size of the 
sample (R. A. Fisher, Statistical Methods, Section 35). Presumably 
your N=n-—8. (ii) The variate x is normally distributed with vari- 
ance o”; hence o?=1/(n—3). (iii) Your formula (1) is identical with 


SN (2-2)? = 3(n—-8) (2-2)? = 3(2-2)?/o?, 


the weighted sum of squares of deviations of z. This quantity is dis- 
tributed as x’, with (k-1) degrees of freedom, under the hypothesis 
that the z’s are all estimates, with the above variances, of a single 
parameter (R. A. Fisher, Statistical Methods, Section 21.03). 

Formula (2) is an approximation to x’, applicable to enumeration 
data. For a discussion, see M. G. Kendall, The Advanced Theory of 
Statistics, Vol. 1, Chapter 12. 

Oscar K&MPTHORNE 
(49) 
QUERY: I am interested in the problem of heredity as an etiological 
factor in homosexuality. 

It was shown by Goldschmidt in his genetic study of butterflies that 
he could establish successive stages of intersexuality from the full male 
through slightly and markedly feminized males to the completely 
transformed female sex intergrade. He also produced a set of inter- 
sexes ranging from complete female through all intermediate stages 
to the male sex intergrade. It cannot be assumed that these results 
apply directly to men but to explain intersexual stages and the trans- 
formation of sex many geneticists have felt that sex is determined by 
sex chromosomes and also by the potentialities of the autosomes. 

It has been shown that changes in temperature, nutrition, light, 
hormones, vitamins, and changes in oxidation rates of anabolism and 
catabolism in organisms (changes in physical and chemical environ- 
ment) may modify the sex of a developing organism to change the 
somal configuration of the organism. 

If the developmental balance theory of Goldschmidt is correct then 
intersexes are all genetically of the female type. If this is true then 
all of the descendants of intersexes would be female. There also 
should be more males found among the siblings of intersexes. If it is 
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true that a percentage of male homosexuals are,real male sex inter- 
grades which are genetically female, but which have lost all morpho- 
logical sex characteristics except their chromosome formula, then an 
undue preponderance of males should be found among the siblings of 
male homosexuals. The normal ratio of 106 males:100 females at 
birth (including live births and stillbirths, excluding stillbirths of 
unknown sex) should then be altered to show an increase in males 
since a series of genotypical females would be among the phenotypical 
male probands. 

I am enclosing the data as determined for two groups. I have at- 
tempted to determine the significance of the difference in the sex ratio, 
but am not satisfied that this method is appropriate to the problem 
or that I am handling it correctly. Does this method determine 
whether or not the differences are significant? What is the standard 
error of the sex ratio? 

Group 1. N=91, males=48, theoretical sex ratio = 106, how rea- 
sonable to get 111? 

p=48/91 = 0.52, q=48 


_ [(0.52) (0.48) 
= 4 i 0.05 


0.52 — 0.50 
0.05 


the o, therefore not significant. 
Group 2. N=148, males = 82, theoretical sex ratio = 106, how rea- 
sonable to get 124? 


p = 106/148 = 0.71, g=1-0.71= 0.29 


o = \/ (148) (0.71) (0.29) = 5.52. 


The difference between 106 and 124 is more than 3 times o, and there- 
fore arises other than by chance. 


= 0.4 times o, which is not 3 times 


ANSWER: In your first group there are two corrections that should 
be made. First, in the formula for o, p = 106/206 = 0.515, the theo- 
retical probability, instead of the sample estimate which you have 
apparently indicated. Second, in the difference between the theo- 
retical and sample probabilities you have put the theoretical value 
equal to 0.5 instead of the correct value shown above. The caleula- 
tions for group 1 are as follows: 


p= 0.515, q = 0.485, = \/(0.515) (0.485) /91 = 0.052 


100 


Sample ratio = 48/91 = 0.527, t = (0.527 — 0.515) /0.052 = 0.28 


In your group 2 the computations seem to have no relation to the 
facts and have led you to a false conclusion. Here, 
p = 106/206=0.515, as before, and «= /(0.515) (0.485) /148 = 0.041 
Sample ratio = 82/148 = 0.554, t= (0.554 — 0.515) /0.041 = 0.95 
As in group 1, the deviation of the sample ratio is not so large as o. 
The sex ratio as you define it has the standard error, 


Loup 
ag l+p n’ 
where p is the hypothetical sex ratio (in your case, p=1.06) and n is 
sample size. 


In groups like yours it may be slightly easier to compute ¢ directly : 
ipa 

Vip’ 

where f is the number of females and m the number of males in the 
sample. 


t 


Grorce W. SNEDECOR 
(50) 

QUERY: A problem has arisen in developing a deer browse sampling 
technique to appraise the available browse supply in deer wintering 
areas in the Adirondack region of New York State and to establish a 
reasonably accurate method of ascertaining the yearly change in the 
browse supply. Because of the extent of the wintering areas, they 
range from 14 to 14 square miles in area, it is not feasible to use random 
samples based on a grid system since the location and appraisal of each 
plot would entail the use of surveying instruments to accurately spot 
a number of scattered plots. 

The method we now employ is as follows: two or three traverse lines 
are drawn at right angles to the main stream bed or valley of each area. 
The lines extend to the perimeters of each area. This insures that all 
altitudinal regions within each area will be represented in approxi- 
mately the same proportion as they exist throughout the entire area. . 
Each traverse line is then divided into 100 equal parts and numbered 
1 to 100. Fifty random numbers are selected from Tippett and the 
numbers drawn are used to spot plots on the traverse line, the plots 
being irregularly spaced due to random selection. The points selected 
are staked and the area within a 25-foot radius comprises the plot to be 
studied. All browse species within each plot are listed on a field 
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schedule and the relative density of each species is noted. Since the 
method employed is a hodge-podge of selective and random sampling 
(selective in choosing traverse lines and random in locating the plots 
along the line) the question arose as to whether it would not be better 
to use selective sampling for both traverse lines and plots along the 
lines. The selection of the traverse line violates the principle of 
random sampling and thus renders the data gained unfit for interpre- 
tative analysis inasmuch as the standard deviation would be meanine- 
less. On the other hand randomization of the plots along the line 
would enable one to calculate the standard deviation and standard 
error but since these statistics would be vitiated by the selection 
originally used in locating the traverse line it is questionable as to 
whether any useful purpose is served in using a random sample of 
plots along the line. 

I will appreciate your views on this subject and will weleome any 
suggestion you may make or any books you may know that treat of 
similar problems. I wish to point out in closing that the main purpose 
of the study is to determine yearly changes in the browse supply. An 
accurate appraisal of the browse supply, or a measurable degree of 
accuracy, on each area is desirable, but probably unattainable without 
random sampling. 


ANSWER: Your two objectives (a) to appraise the available browse 
supply and (b) to ascertain the yearly change in browse supply, indi- 
cate the general field procedure known as the method of double 
sampling, 

For the moment consider objective (a) above. Your present 
method of running traverse lines at right angles to the main stream, 
from ridge to ridge—that is, across the stratification which is roughly 
parallel to the stream—seems efficient and practical from this distance. 
Right here I should like to suggest that the number of sample points 
along the traverses (centers of sample-plot circles of 25-ft. radii) be 
proportional to length of traverse—such, for instance, as one in each 
100 feet of traverse length, taken systematically. 

If you are required to appraise the browse supply of each wintering 
area with satisfactory precision, I imagine that 2 or 3 traverse lines as 
now employed are insufficient. Perhaps a random sample of 2 tray- 
erses to each mile of distance, running at right angles to an imaginary 
line from the down-stream end to the top of the valley, would prove 
workable. This scheme implies the division of a wintering area into 
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blocks having one dimension of one mile, but the other dimension vary- 
ing with distance from base line to ridge. It is discussed in detail in 
Chapter VII of ‘‘Sampling Methods in Forestry and Range Manage- 
ment’’ by Schumacher and Chapman (Bulletin 7, School of Forestry, 
Duke University, 1942). 

On the other hand, if you are required to appraise the supply of 
only the aggregate of all areas, and not of each individual area, the 
field work may be restricted considerably, perhaps to your present 
scheme of 2 or 3 traverses in each wintering area. On the supposition 
that you require the sampling error of your appraisal, a random sample 
of traverses should be taken from each area. ‘ 

The above discussion has to do with objective (a) only—that is, the 
appraisal of the browse supply of the particular year. The adjust- 
ment to the supply of the succeeding year should require relatively 
little field work. For this purpose, the circular plot of 25-ft. radius 
should serve as a satisfactory sampling unit. If in the first year a 
certain (relatively small) number of these are identified by a perma- 
nent post at their centers, and in the second (and each succeeding) 
year, the browse supply on these permanent plots alone be measured, 
then the regression of the second-year browse on the first-year browse 
of these plots is the adjusting equation required. One inserts into the 
equation the average browse supply (on the circular plot basis) as 
obtained in the comparatively heavy sampling of the first year and 
solves for the corresponding value of the second year. 

It thus seems that a somewhat heavy sampling need be made only 
once every several years—that is, only at such time as the regression of 
succeeding years’ measurements of individual plots on those of the 
base year no longer serve as satisfactory adjustment. 

F. X. SCHUMACHER 
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ABSTRACTS* 


(40) 


GREENWOOD, MARY L. (University of Connecticut). Studies 
on the Palatability of Kale in Relation to Variety and Cooking Pro- 
cedure, 

The relative palatability was tested by means of consumer prefer- 
ences. The ranks were converted to mean normal deviates. These 
transformed scores were then compared in a series of analyses of 
variance. Where a larger number of alternatives had to be compared 
than could be handled conveniently at one time, the experiments were 
arranged in balanced incomplete blocks. The method of ranks seems 
adapted to relatively complex experiments for measuring food prefer- 
ences. 


(41) 
BROWN, GEORGE W. (Iowa State College). Discriminant 
Functions. 

The logical foundations of discriminant function techniques are 
developed from an elementary point of view. The basic model postu- 
lates underlying ‘‘measurements’’ of the form 

gy Ea (6 bap Bw + E34 


E=aet€e 


Ep = Op + Ep 
where the a; and 8 are unknown constants, ¢; are independent errors 
of measurement, and w is a parameter varying over a set of popula- 
tions. In addition, it is postulated that the observables are of the form 


Dp 
Cm Solan en Map lay Dp 
n=1 


where the {Im,} are unknown coefficients. If {omn} is the matrix of 
variances and covariances of the x’s, and {o™"} is the inverse of this 
matrix, and if b, is the regression coefficient of am on w, then the 

* All papers read at the joint meetings of the Biometrics Section and the Institute 


of Mathematical Statistics held in conjunction with the 113th Annual Meeting of the 
American Association for the Advancement of Science in Boston, December 27-29, 1946. 
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Dp 
function X=)) o”"bnt, is proportional to the hypothetical mea- 


mna1 

surement £,, the only one of the &n which is related to w. This X dis- 
eriminates among w values to the same extent that £,, were it available, 
would discriminate. X is, of course, also that linear combination 
(unique except for a constant of proportionality) of «’s for which the 
ratio of the mean square for regression to the mean square for error 
is a maximum, or what is the same thing, X is the linear combination 
of x’s which has greatest correlation with w. The coefficients of X are 
proportional to those which would have been obtained from a formal 
regression analysis of w on 21, %2,..., Up The usual procedure for 
testing against subsets of the x’s still holds in this situation, in spite of 
the inversion of roles of dependent and independent variables. 

Two generalizations to multi-dimensional population parameters 


are possible. Let wi, We,.--, Ws, be parameters associated with a 
set of populations. The case s=2 is considered. Let there be hypo- 
thetical measurements &1, &, . - - &p, with 


Ey = 01+ Biwt yi + &1 
E> = Oe + Bol + YoU + €2 
&3 = 3+ €3 


Ep = Opt Ep 
Let b», be the regression coefficient of %m on U, Cm the regression coeffi- 
cient of m on v, and let omn be the variance-covariance matrix among 
the {em}. Then X1=Lo"™bmr, and X,=DXo™Cn&m are independent 
linear combinations of &, and &, alone, so that x1 and X, together con- 
tain all the information in the «’s relevant to wv and v. X, and Xz 
together yield the best 2 -dimensional description of the populations. 
The other generalization corresponds to finding the best 1-dimensional 
description of the populations, leading to the solution of a determinan- 
tal equation. X, and X, are the 9-dimensional discriminant function ; 
the best 1-dimensional linear function, in the sense of ratio of mean 
square for regression to mean square for error, is the ‘‘reduced’’ 
discriminant function. In general, when s > 2, there is a hierarchy 
of reduced discriminant functions with dimensionality 1, 2,...,s—1. 
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NEWS AND NOTES 


Statisticians and those interested in statistics from the United 
States and foreign countries will convene at Virginia Polytechnic 
Institute, Blacksburg, Virginia, for the Statistical Summer Session, 
August 5 to September 5. This Summer Statistical Session is being 
sponsored jointly by Virginia Polytechnic Institute, University of 
North Carolina, University of Michigan, Iowa State College, and the 
Federal Bureau of Agricultural Economics, The faculty will con- 
sist of: WALTER A. HENDRICKS, B.A.K., U.S.D.A.; RENSIS 
LIKERT, University of Michigan; H. L. LUCAS, University of 
North Carolina; MAURICE G. KENDALL, England; GEORGE 
W. SNEDECOR, Iowa State College; FRANK YATES, Rothamsted 
Experiment Station, England; EARL E. HOUSEMAN, B.A.E., 
U.S.D.A.; RAYMOND. J. JESSEN, Iowa State College; and BOYD 
HARSHBARGER, Virginia Polytechnic Institute. The following 
courses will be offered for credit: Engineering Statistics; Statistical 
Methods; Design of Animal Experiments ; Schedule Design and Inter- 
view Techniques for Sample Surveys ; Sampling Design and Analysis ; 
Mathematical Theory of*Sampling; Seminar ; Mathematical Statistics; 
and Experimental Design. In addition to the faculty, probable 
Seminar speakers are: W. F. CALLENDAR, W. G. COCHRAN, 
GERTRUDE M. COX, W. E. DEMING, GEORGE GALLUP, MOR- 
RIS HANSEN, HAROLD HOTELLING, ARNOLD KING, and 
CHARLES F. SARLE. Inquiries regarding the Summer Session 
should be addressed to Boyd Harshbarger, Professor of Statistics, 
Statistical Summer Session, Virginia Polytechnic Institute, Blacks- 
burg, Virginia. . . . The Biological Methods Group of the Society of 
Public Analysts and Other Analytical Chemists held its annual meet- 
ing on December 16, 1946. Two papers were read: 1. ‘‘The assay of 
anti-thyroid substances using tadpoles”’ by H. M. BRUCE. In this 
paper she was concerned almost entirely with matters of biological 
technique. 2. ‘‘The computation of microbiological assays of amino- 
acids and other growth factors’? by ERIC B. WOOD. He points out 
that microbiological assays which did not show linear relationship 
between dose and response are nearly always linearized when the 
logarithms of both the dose and the response are taken as metameters. 
..-R. HOWARD PORTER, Director, Instituto Agronémico Na- 
cional, left Iowa State College a year ago for Paraguay to work on 
the Food Production Program. His work is primarily concerned with 
varietal testing of many kinds of crops and seed production of the 
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improved varieties which have been tested in the past. In addition, 
he is doing work with seed disinfectants and agronomic experiments 
such as spacing and fertility tests. There is a nursery for testing 
many kinds of vegetables for their adaptation to this region... . 
HOWARD L. STIER has returned to the University of Maryland as 
Professor of Marketing and head of the department after more than 
five years in the army. He had been a member of the Department 
of Horticulture at the University before entering the army in 1941. 
During the war, Dr. Stier served as Chief Statistician of the Quarter- 
master Corps and later as Director of the Statistics and Progress 
Reporting Division of War Assets Administration, rising from the 
rank of Ist Lieutenant to Colonel. He was awarded the Legion of 
Merit for ‘‘outstanding service’? in the Quartermaster Corps... . 
BURCH H. SCHNEIDER, Professor of Animal Husbandry, West 
Virginia University, Morgantown, deserves a word of commendation 
for his unlimited patience in the compilation of the material for ‘‘ Feeds 
of the World.”?. ..T. A. BANCROFT, now at the University of 
Georgia, Athens, will go to Alabama Polytechnic Institute at Auburn 
the middle of June. He goes to organize and direct a new statistical 
laboratory under the sponsorship of the Dean of the Graduate School 
and the Alabama Polytechnic Institute Research Council. The func- 
tions of the Statistical Laboratory are: (1) to provide instruction in 
statistical methods on the senior and graduate level, (2) to provide 
a consultation service for research scientists at the university and in 
the state, and (3) to provide computational service for college and 
state agencies. Both the State Agricultural and Engineering Experi- 
ment Stations are located at Auburn. ALONZO COHEN, now work- 
ing in mathematical statistics at the University of Michigan, will take 
his place at the University of Georgia. . . . MARGARET E. RICE 
has resigned her position as statistician with the TVA at Wilson Dam 
to become supervisor of public health statistics at Mississippi State 
Board of Health at Jackson. She assumed her new position February 


1.... The Tennessee Public Health Association held its annual 
meeting in Nashville, May 5-7. ANN DILLON arranged the pro- 
gram. ... W. EDWARDS DEMING, the man who travels to Eng- 


land, India, and Japan, should provide some highlight reports for us, 
especially when he writes, ‘‘I had a wonderful time in India and came 
away with many puzzling thoughts.”’ . . . E. J. WILLIAMS, Division 
of Forest Products, Council for Scientific and Industrial Research, 
South Melbourne, Australia, writes: ‘‘Since I began work here at the 
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Division of Forest Products some six years ago, the statistical work 
required by the Division has increased considerably. There are now 
two of us on the statistical staff, with five computers. Much of the 
work done here is routine, such as computing of test results, and 
analyses of variance. However, we get a number of interesting prob- 
lems. Timber being subject to defect and loss, the analysis of experi- 
ments involves quite a lot of treatments of non-orthogonal data, for 
which we have been developing some useful approximate methods, as 
well as iterative methods for arriving at the exact solution. We have 
come across a number of problems in design, such as factorial experi- 
ments, which have to be fitted into unbalanced incomplete block de- 
signs, so as to confound only interactions. On the theoretical side, I 
have been playing round with significance tests in multivariate analy- 
sis, but so far do not see any exact approach. R. A. FISHER, I be- 
lieve, has been delving into the problem. You will be interested to 
hear we are making use of 7x 7 Latin square designs in studying the 
effect of two consecutive processes in laboratory sheet making ; time 
of drawing of stock, and position in the press. By randomizing the 
sheets in such a way that each sheet from one time of drawing occurred 
at a different one of the seven positions in the press, in successive repli- 
cations, we were able to disentangle the two effects satisfactorily.”? He 
inquired whether we have met David Duncan. We have. The Edi- 
torial Committee would be glad to have others send us brief reports 
regarding their work. ...H. J. MILLER, Assistant Professor of 
Plant Pathology at Pennsylvania State College, comments on this bul- 
letin : ‘‘For the average biological worker, it affords an opportunity to 
have questions answered by competent statisticians such as are not 
available in many of our institutions.”’? . . . C. H. GOULDEN, Officer 
in Charge, Cereal Division, Department of Agriculture, Winnipeg, 
Canada, writes: ‘‘I should like to take this opportunity of stating that 
the Biometrics Bulletin is filling a very definite need. It is not pos- 
sible for most of us to read more than one half of the papers that are 
being published at the present time on statistics. There must be some 
journal that either digests articles from other journals or presents gen- 
eral articles by well-known writers on particular phases of statistical 
methods.’’ Mr. Goulden thinks our new format is a real improve- 
ment. We hope you agree. 
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